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ABSTRACT 


investigation was performed to determine the feasibility 
of using a two plate model to approximate the Tracking and Data 
Relay Satellite (TDRS) in orbit propagation, taking into account 
the effects of solar radiation pressure. The two-plate model 
comprises one plate which always points to the earth, and the 
other which is hinged to an axis normal to the orbital pland and 
is always rotated so that its normal makes a minimum angle wit 
the direction of the sun. 

The results indicate that it is sufficient to take three par- 
ameters (i.e. , the areas of the two plates and the reflectivity 
of the earth-pointing plate) to achieve an accuracy of one meter 
during a 24-hout orbit propagation. 
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SECTION 1 ~ INTRODUCTION 


Most of the work involving solar radiation pressure on 
orbiting satellites has so far been limited to those which are 
spherical and which have circular nominal orbits. The compara- 
tively few studies which are less restrictive are still based 
on very simplified models. For example, the work of Eliasberg^^ 
deals with an elliptically orbiting spherical satellite and is 

concerned with the first order perturbation effects expressed 

1 2 ) 

m terms of Keplerian elements. The work of Fang deals 
with a circularly orbiting spherical satellite with a perfectly 
reflecting earth-pointing disk. It is concerned with the first 
order effects expressed in terms of along-track, cross-track 
and radial components. Moreover, it also deals with the physical 
insights into the modeling errors connected with tracking and orbit 
determination of the Tracking and Data Relay Satellites (TDRS) . 

On the other hand, the work of Georgevic ' deals only with the 
computation of solar radiation force on a cylinder and on a 
parabolic reflector, but does not deal with an orbiting satellite. 
(Even then, the computation for the parabolic reflector is further 
simplified by assuming that the ratio of the force on the illumin- 
ated area of the reflector and the force on the whole area of the 
reflector is the same as the ratio of the corresponding projected 
areas. It is obvious that this assumption, introduced to eliminate 
the cumbersome self-shadowing effects, is not really correct.) 

The present work is concerned with the solar radiation 
effects on the TDRS illustrated in Figure 1.1, and modeled 
as comprised of the 69 components listed in Table 1,1. (In 
the course of the present study, a novel method, simple in 
comparison to other existing methods, for computing self-shadowing 
was formulated but this consideration was not included in computing 
the net solar radiation force on the satellite) . The orbit of 
the TDRS is taken to be representative of a realistic one in that 
it is not nominally perfectly circular. The study also considers 
the question of how accurately the 69-component TDRS can be 
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approximated by a two-plate model, i.e., one plate is hinged to 
an axis normal to the orbital plane and is always rotated to 
make a maximum angle with the sun, while the other plate is 
always earth-pointing. This two-plate model has the capability 
of handling up to four solve-for parameters, i.e., the area and 
reflectivity of each of the two plates. 

Section 2 is concerned with the analysis of a differential 
correction procedure to obtain the values of these four solve-for 
parameters. A reference orbit for the 69-component TDRS is 
first generated. Its orbital position at regular intervals of 
one hour is then used as epoch elements of the two-plate model 
to obtain the values of the parameters which yield the best 
approximating orbit over the next 24 hours. 

Section 3 summarizes the numerical results obtained in 
this feasibility study, and presents tabulated and graphical 
results for rapid comparisons. 

Section 4 discusses the quality of the results, and the 
applicability of the two-plate model for use in orbit determination 
purposes . 
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TABLE 1 . 1 



DETAILS OF 

69 -COMPONENT MODEL 


Component v 

Description 

Area (m^) 

Reflectivity 

Components of Normal 
X Y Z 

1 

Solar Panel 1 

14.7553 

0.0 

Normal makes 
minimum angle 
with sun 

2 

Solar Panel 2 

14.7553 

0.0 

ti «• «* 

3 

SGL Antenna 

3.13761 

0.5 

Points to 
White Sands 

4 

C-Band Antenna 

2.67112 

0.5 

Points to 
Los Angeles 

5 

Solar Sail 

0.90593 

1.0 

0 0 1 

C 

Antenna Feed 
1 (top) 

0.00462 

1.0 

0 0 1 

1 

Antenna Feed 
2 (top) 

0.00462 

1.0 

0 0 1 

t 

Antenna Feed 
1 (bottom) 

0.29570 

1.0 

0 0-1 

9 

Antenna Feed 
2 (bottom) 

0.29570 

1.0 

0 0-1 

10 

Antenna Feed 
1 (side) 

1.11771 

1.0 


11 

Antenna Feed 
2 (side) 

1.11771 

1.0 


12 

Main Body (top) 

4.03665 

0.5 

0 0 1 

13 

Main Body 
(bottom) 

4.03665 

0.5 

0 0-1 

14 

Main Body 
(side 1) 

1.1599 

0.5 

0 10 

15 

Main Body 
(side 2) 

1.1599 

0.5 

.866 .5 0 

16 

Main Body 
(side 3) 

1.1599 

0.5 

.866 -.5 0 

17 

Main Body 
(side 4) 

1.1599 

0.5 

0-10 

18 

Main Body 
(side 5) 

1.1599 

0.5 

-.866 -.5 0 

19 

Main Body 
(side 6) 

1.1599 

0.5 

-.866 .5 0 

20-69 

Sections of 

Commuted by 

1.0 

Computed by POLYK 


stationary 

antennae 
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SECTION 2 - ANALYSIS 


This section is concerned with the analysis of a 
differential correction procedure to obtain the values of 
the four solve-for parameters incorporated into the two-plate 
model . 


For convenience, let us introduce the following 
notation: 


f (t) = position vector of 69-component TDRS 
at time t 

g(t,5<) = position vector of two-plate model at time t 
1* s (a*, , tt’*,**'*) = 4 parameters for two-plate model 

2-5 4 ^ . j 

Q = loss function defined as Z Htt)- 3fc,s)| 

L-l 


The problem is then to obtain the values of « such that Q 
is a minimum. It is obvious that the minimum of Q is given by 
the necessary condition = 0 where j = 1, 2, 3, 4. 

The loss function Q may also be written as 

25 


a - 1 [fto-jte,*)]- [foo- at*,*)] 

Z*i 


( 2 - 1 ) 


The necessary condition for minimum is 


• fa '" ) = 0 3 = 1 ' 2 ' 3 ' 4 

J C»I J 


( 2 - 2 ) 
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The function g^t^c*) may be expanded in a Taylor series about 

—A 

an a priori value <X t . 


t Ck,*) - + i ^ 


A« k + 


k* I 


(2.3) 


0< - Pit 


where Ao<, is defined as 
k 


Ao< v = c* k - <x c ^ k 


k=l ,2,3,4 (2.4) 


Let 


(2.5) 


Substitution of equations (2.3) and (2.5) into (2.2) yields 
to first order 

z {!>*&)- cl ^«oi • } = ° 


l-l 


Ic-, 


“J 1 o(= « c 

j = l , 2 , 3 , 4 (2.6) 


Interchanging the order of summation yields 

4 - ze - »?(*:,*) I . ■Safe.sj 


z(z[ 

k*I t=l 




3« 


UK 


k 


25 r- 

- z 

£*» 


tgCt:,*) 


3* 


£:,«.)] 


J ** 


Let a. 


y r [ # (t ^g) 1 

jk - Z L a«; im ~ 9 « k imt >J 


(2.7) 


j=l, 2,3,4 


( 2 . 8 ) 


d — j r k- 1 / 2 f 3 t 4 


i»j 


25 


= ZL 

l=M 


r *) 


■ 


I <** <*1 
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(2.9) 


Equation (2.7) then becomes 

4 

2 *J k ~ h j = l / 2 , 3 , 4 (2.10) 

l<-| 

It remains to solve for where k=l y 2 f 3,4. This constitutes 

the first iteration in the differential correction procedure to 
solve for the value of o( such that Q is a minimum. 
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SECTION 3 


RESULTS 


This section summarizes the numerical results in this 
feasibility study, and presents tabulated and graphical results 
for rapid comparisons. 

Numerous computer runs were made for the TDRS with epoch 
elements : 

X c = 31,662,513.0 m 

Y q = -27,523,890.0 m 

Zq — 0.0m 


X Q = 2,012.15997 m/s ec 

Y q = 2,314.7253 m/s ec 

Z = 376.58528 m/sec 
o 

One set of runs had epoch time set at Day 183.0, Year 1980 (i.e., 
July 1, 1980 which is close to the summer solstice) , and another 
set had epoch time set at Day 275.0, Year 1980 (i.e., October 1, 
1980 which is close to the autumnal equinox) . Each of these sets 
of runs was made for the cases of N=l,2,3 and 4 parameters. The 
parameters were aligned in the following sequence which is 
probably the order of decreasing importance: 

2 

a = area of sun-pointing plate (M ) 

2 

= area of earth-pointing plate (M ) 

= reflectivity of earth-point plate 
= reflectivity of sun-pointing plate 

The initial values of these parameters used in the first iteration 
of the differential correction procedure were taken to be the 
following : 


N 

= 1 

2 

3 

4 

P 

o 

= 36.6 

29.5 

29.5 

29.5 

a 

0,2 

= 0.0 

18.81 

18.81 

18.81 

ot _ 

= 0.0 

0.0 

0.74 

0.74 

0,3 

a 0,4 

= 0.0 

0.0 

0.0 

0.0 


The tolerance e for testing convergence of the iterations 
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(i.e. , \&tt |<€) was taken to be 0.01. The values for the stepsize 6 
for computing the partial derivatives was taken to be 0.1. In the 
runs described above, convergence was achieved after two or three 
iterations. Figures 3.1 and 3.2 summarize the results of these 
computer runs. The symbols appearing on the horizontal axis in 
these figures have the following connotations : 

0 = The sun is overhead with respect to the satellite 
U = The sun is underfoot with respect to the satellite 
A = The satellite is moving away from the sun 
T = The satellite is moving toward the sun 
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Fiquro 
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Figure 


SECTION 4 - CONCLUSION 


unis section is concerned with a brief discussion of the 
results obtained in Section 3, and also an attempt to relate 
them to results obtained in other investigations. 

From Figures 3.1 and 3.2, it is seen that the following 
observations may be made : 

1. The maximum deviation of the two-plate model from 
the 69-component TDRS is essentially cyclical with 
a 12-hour period. 

2 . The magnitude of these deviations decreases as the 
number N of parameters increases , as would be 
expected . 

3. There is a pronounced phase change of this period 
curve in going from the case of N=1 to N=2. 

4 . The curves for the cases of N=3 and N = 4 are essentially 
identical. Moreover, the amplitude of oscillation is 
so small that they are almost constant. 

A computer run was also made comparing the 6 9 -component TDRS 
with and without solar radiation pressure effects. The results 
are plotted in Figure 4.1. 


It is noted that in this case the curve is essentially sinusoidal, 
unlike those in Figure 3.1 and 3.2. Moreover, the maximum deviation 
occurs when the sun is overhead or underfoot in Figure 4.1, unlike 
the cases of N=1 and N=2 in Figures 3.1 and 3.2. 

Finally, it is interesting to recall the results obtained by 

( 2 ) 

Fang ' ' who observed that : 

1. In orbit propagation, the least perturbation occurs 
when the sun vector is parallel to the satellite 
velocity vector in the beginning , and the worst 
perturbation occurs when the sun is overhead or 
underfoot in the beginning . 

2. TDRS orbits determined from a one-day tracking arc 
tend to be less sensitive to solar pressure errors 
if the tracking arc begins when the sun is directly 
overhead or underfoot. This is contrary to the 
(previous) result for solar pressure perturbations 
in the absence of tracking. 


1-13 



1-14 


Figure 


From the numerical results of this feasibility study, it is 
seen that the 69-component TDRS can be accurately replaced by 
the two-plate model. It suffices to take only three parameters 
to achieve an accuracy to within about one meter. Moreover, it 
is sufficient to use only one approximating orbit throughout 
the 24 hour period, instead of 24 approximating orbits regularly 
spaced throughout the day . 


REFERENCES 


(1) Eliasberg, P.E., "Introduction to the Theory of Flight of 
Artificial Earth Satellites," Moscow (1965). Translation, 

— n / i n r -7 \ 

jtlUbaicm (i^u / j • 

(2) Fang, B.T., "Short-Term Solar Pressure Effect and GM 
Uncertainty on TDRS Orbital Accuracy A Study of the Inter 
action of Modeling Error with Tracking and Orbit Determina- 
tion," EG&G/Washington Analytical Services Center, Inc., 
Report No. 002-78 (1978) . 

(3) Georgevic, R.M. , "The Solar Pressure Force and Torque Model, 

J. Astron. Sc., Vol. XX, No. 5, pp. 257-274. (1973). 


1-15 


